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ABSTRACT 



Context. Discovery of the 6.7-hour periodicity in the X-ray source IE 161348-5055 in RCW 103 has led to investigations of the 
nature of this periodicity. 

Aims. To explore a model for IE 161348-5055 wherein a fast-spinning neutron star with a magnetic field ~ 10 12 G in a young pre- 
Low-Mass X-ray Binary (pre-LMXB) with an eccentric orbit of period 6.7 fir operates in the "propeller" phase. 
Methods. The 6.7-hr light curve of IE 161348-5055 is modeled in terms of orbitally-modulated mass transfer through a viscous 
accretion disk and subsequent propeller emission. Formation of eccentric binaries in supernovae and their subsequent tidal evolution 
are studied. 

Results. The light curve of IE 161348-5055 can be quantitatively accounted for by models of propeller torques of both Illarionov- 
Sunyaev type and Romanova-Lovelace et al. type, and spectral and other properties are also in agreement. Formation and evolution 
of model systems are shown to be in accordance both with standard theories and with X-ray observations of IE 161348-5055. 
Conclusions. The pre-LMXB model for IE 161348-5055 and similar sources agrees with observation. Distinguishing features be- 
tween this model and the recently-proposed magnetar model need to be explored. 
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1. Introduction 

The point soft X-ray source IE 161348-5055 (henceforth IE) 
near the center of the young (~ 2000 yr old) supernova remnant 
(SNR) RCW 103 has attracted much attention lately, following 
the discovery of a strong 6.67 hr periodic modulation in IE by de 
Luca et al. (2006, henceforth dL06) from a deep XMM-Newton 
observation of the source in 2005. IE was discovered in 1980 
( |Touhy & Garmire 1980[ ) as a soft Einstein X-ray source. The 
original interpretation as an isolated neutron star was found to 
be untenable in view of subsequent discovery by other X-ray 
satellites {e.g., ROSAT, ASCA, Chandra) of the large variability 
of IE on the timesacle of a few years (dL06). A periodicity at 
~ 6 hr was first hinted at by Chandra observations, but the first 
clear, strong detection came from the above 2005 observations 
of dL06, who also showed the existence of this periodicity in the 
data from earlier 2001 observations of IE with XMM-Newton, 
when the source luminosity was higher by a factor ~ 6 during the 
course of its sequence of several-year timescale outbursts nen- 
tioned above, documented by these authors from archival data. 

The nature of the above 6.67 hr periodicity is an interest- 
ing question, on which preliminary discussions were reported 
in dL06. Recently, Pizzolato et al. (2008, henceforth P08) have 
proposed a model for the IE system wherein it is a close binary 
consisting of a magnetar, i.e., a neutron star with a superstrong 
magnetic field ~ 10 15 G, and a low-mass companion. The 6.67 
hr periodicity is identified in this model with the spin period of 
the neutron star, to which this young neutron star has been spun 
down in such a short time by the torques associated with its enor- 
mous magnetic field. This period has also been proposed by P08 



to be in close synchronism with with the orbital period of the bi- 
nary, in analogy with what is believed to be happening in Polar 
Cataclysmic Variables or AM Her-type systems. The observed 
X-ray emission from IE is that from the magnetar in this model. 

In this paper, we explore an alternative model for the IE 
system wherein it is a close binary system consisting of a young 
neutron star with a canonical magnetic field ~ 10 12 G, and 
a low-mass companion, i.e., a pre-Low-mass X-ray Binary 
(henceforth pre-LMXB), such as are believed to be the standard 
progenitors of Low-mass X-ray Binaries (henceforth LMXBs). 
Such pre-LMXBs are born after the common-envelope (CE) 
evolution phase of the original progenitor binary system con- 
sisting of a massive star and a low-mass companion, which 
leads to the formation of a binary consisting of the He-core 
of the original massive star and the low-mass companion 
dGhosh 20071 and references therein). The He-star susequently 
explodes in a supernova, leading to a neutron star in orbit 
with a low-mass companion, i.e., the pre-LMXB referred to 
above. This is the standard He-star supernova scenario for the 
formation of LMXBs ( Gh osh 20071 and references therein). 
The 6.67 hr periodicity is identified in our model with the 
orbital period of the binary. In our model the young neutron 
star is still spinning very rapidly, with a canonical spin period 
~ 10 - 100 ms, and is operating in the "propeller" regime, 
wherein any matter approaching the fast-rotating magne- 
tosphere of the neutron star is expelled by the energy and 
angular momentum deposited into it through its interaction with 
the magnetospheric boundary (Illarionov & Sunyaev 1975 
henceforth IS75, IDavies etal. 19791 



|Illarionov & Kompaneets 1990] 
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Mineshige et al. 1991 
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llllarionov et al. 19931 [Ghosh 1993] and references therein, 
henceforth G95, [Lovelace et al. 19991 hence forth LRB99 
IRoma nova et al. 2004] IRomanova et al. 20051 henceforth 



RUKL05, Ustyugova et al. 2006| henceforth UKRL06). 



The observed X-ray emission from IE in our model is that 
from the propeller: indeed, it is well-known that soft X-ray 
transients (SXRTs) like Aquila X-l and others (see Sec. 17. It 
go through low/quiescent states during the decay of their out- 
bursts, during which their luminosities and spectral properties 
are very similar to those of IE, and the neutron stars in them 
are believed to be in the propeller regime (Campana et al. 1998 
Stell a et al. 20001 1. The observed 6.67 hr periodicity in our model 
is due to the orbital modulation of the supersonic propeller, 
which is caused by the orbital modulation of the mass-transfer 
rate in the eccentric binary orbit of a young system like IE. It 
is well-known that young post-SN binaries with low-mass com- 
panions like IE are almost certain to have eccentric orbits, due 
to the large eccentricities produced in such systems in the SN 
explosion (see Sec. 16. It and the duration of the subsequent tidal 
circularization compared to the ages of systems like IE (see Sec. 
I6.21 i. By contrast, SXRTs are believed to be old LMXB systems 
with circular orbits, where such modulation will not occur. 

We show in this work that the 6.67 hr light curve of IE can be 
accounted for quantitatively by our model for propeller torques 
of both Illarionov-Sunyaev type and Romanova-Lovelace et 
al. type (see Sec. and that the observed spectral and other 
characteristics are also in general agreement with our overall 
picture. Thus, further diagnostic features need to be explored in 
order to distinguish between our model and the magnetar model 
as a viable description of this and similar sources. 

2. Propeller phase in pre-LMXBs 

In a pre-Low-Mass X-ray Binary (pre-LMXB: see above), 
the newborn, fast-rotating neutron star is unable at first to 
accrete the matter that is being transferred from the com- 
panion through the inner Lagrangian point L\, because of 



the fast rotation of the neutron star (IS75, IDavies et al. 1979 
[T98T1 |Illarionov & Kompaneets 1990| I Mineshige et al. 1991 
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Illari onovefq/. 19931 G95, LRB99, [Romanova et al. 2004 



RUKL05, UKRL06). Because of its large angular momentum, 
this matter forms an accretion disk and reaches the magneto- 
spheric boundary of the magnetized neutron star, whereupon this 
ionized matter interacts with the fast-rotating neutron star's mag- 
netic field, and the energy and angular momentum deposited in 
it by magnetic stresses associated with this fast-rotating mag- 
netic field expel it. This is the propeller phase of the system 
(IS75), during which the neutron star spins down as it loses an- 
gular momentum and rotational energy. During this propeller 
phase, the disk matter at the magnetospheric boundary is shock- 
heated as the "vanes" of the supersonic propeller (IS75) hit it, 
and the hot matter emits in the soft X-ray band. This emission 
appears unmodulated at the neutron-star spin frequency (as op- 
posed to the X-ray emission from canonical accretion-powered 
pulsars, which comes from the neutron-star surface) to a dis- 
tant observer, who sees only the total emission from the heated 
matter at the magnetospheric boundary. Observations of tran- 
sient low-mass X-ray binaries {i.e., the soft X-ray transients or 



SXR Ts) like Aquila X-l ( |Campana et al. 1998) and SXJ 1808.4- 
3658 (Stella et al. 2000) in quiescence, when the neutron stars in 
them are thought to be operating in the propeller phase, amply 
confirm this point. 

The propeller luminosity L during the above phase is given 
by L — Nat, where N is the propeller torque acting on the neutron 



star and u> is its spin angular velocity. The propeller torque N 
was first estimated by IS75 in their pioneering suggestion of this 
mechanism, and subsequent work over approximately the next 
two decades considered variations of this torque under different 
circumstances, as summarized in G95. These works addressed 
themselves largely to quasi-spherical accretion, and we shall call 
this kind of propeller torque the Illarionov-Sunyaev type (or IS- 
type for short) torque, which was widely used in that time-frame 
in propeller spindown calculations. In the 2000s, Romanova, 
Lovelace and co-authors reported a series of calculations of the 
propeller effect for disk-accreting magnetic stars based on their 
numerical MHD simulations (Romanova et al. 2004, RUKL05, 
UKRL06; also see the analytic estimates in LRB99). We shall 
call the propeller torque obtained from this line of work the 
Romanova-Lovelace et al. type (or RUKL-type for short) torque. 
In this work, we shall consider both IS -type and RUKL-type pro- 
peller torques for the problem at hand. 

Consider IS-type torques first. For such fast-rotating neutron 
stars as we are concerned with in this work, the propeller oper- 
ates in the supersonic regime, and its torque is given by (G95 
and the references therein), 



_ 1 pW Q(r m ) 
6 GM X o) 



(1) 



In Eq. (HJ, Q(r,„) is the Keplerian angular velocity at the mag- 
netospheric radius r m , p is the magnetic moment of neutron star 
and M x is its mass. Combining this equation with the standard 
expression for the magnetospheric radius (IGhosh 2007K viz., 



M V2GM.v 



(2) 



where M is the rate at which transferred matter arrives at the 
magnetospheric boundary, we obtain the following expression 
for the propeller luminosity: 



L35 ~ 5Mj 4 (P spin /0.ls)- 2 ixl m; 



(3) 



In Eq. ([3), Mi 4 is M in units of 10 14 g s L35 is L in the units 
of 10 35 erg s _1 , P S pin 1S the neutron-star spin period, /Z30 is the 
neutron-star magnetic moment in units of 10 3() G cm 3 , and m x 
is the neutron-star mass in units of solar mass. As neutron stars 
are thought to have P sp i„ ~ 0.01 - 0.1 s at birth, and as the 
propeller phase is thought to end when the spin period is longer 
than P sp i„ ~ 0.1 - 1 s, we have made the canonical choice for the 
expected scale of P S pm m systems like IE. Equation (0 clearly 
shows how the propeller luminosity scales with the mass-arrival 
rate M, and essential neutron-star properties, namely, its spin 
period P sp m, its magnetic moment p, and its mass M x . 

Now consider RUKL-type torques. These authors summa- 
rized the results of some of their extensive MHD simulations in 
RUKL05 and UKRL06 in terms of power-law fits to these re- 
sults, showing that the scaling of the total propeller torque N 
with the magnetic moment p and the spin rate a> of the neutron 
star was 



N oc p lA OJ 2 . 



(4) 



However, the scaling of N with M was not available from the 
above references, because only the parameters p and u (and also 
the turbulence and magnetic diffusivity parameters of the disk: 
see below) seem to have been varied in the series of simulations 
reported in these references. In order to estimate the scaling of N 
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with M for RUKL-type torques, we proceeded in the following 
way. 

First, we did an analytic estimate in the following manner. 
In their analytic study, LRB99 argued that the radius of the 
inner edge of the disk should depend on the stellar rotation rate 
to in addition to the parameters p and M that r m (see above) 
depended upon. The scaling with to, p, and M that these authors 
derived was revised in UKRL06, the final result being given as 
oc p}l 2 M~ l ^ur' 1 ^. (Note the closeness of the scalings with p 
and M with those which apply to r m , as given above.) 

In a simple first approach, if we argue that a reasonable es- 
timate of the torque scalings may be obtained by replacing r m 
with in Eq. ([TJ for disk accretion, we arrive at the scaling 



N CC p 



S/4,.,11/8^3/8 



(5) 



for RUKL-type torques. Noticing the qualitative similarity of the 
the scalings with p and to in Eq. (O with those of the actual 
RUKL-type torque given in Eq. |@}, and furthermore the quan- 
titative closeness for the scaling with p, we argued that the best 
estimate would be to use the scalings of Eq. © for p and a>, and 
the scaling of Eq. © for M, thus arriving at a suggested scaling 
for the RUKL-type torque as 



Noc/j 1 



l co 2 M 3/s . 



(6) 



Before proceeding further, we recognized that RUKL-type 
torques may arise from more complicated interactions than are 
describable by the above arguments, and so attempted to ver- 
ify the above M scaling by further comparison with RUKL re- 
sults. To this end, we noted the correlated variations of N and M 
recorded in Figure 4 of Rom ano va et al. 20041 and fitted the two 
prominent peaks in N and M at the extreme right of this figure 
to a power law. This gave an exponent m 0.37, coincident with 
that in Eq. (|6]i within errors of determination. With this support, 
we use the scalings of Eq. (O for the RUKL-type torque in this 
work, deferring further considerations to future publications. 

In order to obtain the dimensional values of the RUKL-type 
propeller torques and related varaiables, we now insert the ref- 
erence units for the RUKL simulations given in RUKL05 and 
UKRL06, thus obtaining for the torque: 



N33 



14 



(7) 



Here, Afo is the propeller torque in units of 10 g cm s , the 
units of other variables are as before, and we have kept the values 
of the turbulence and magnetic diffusivity parameters of the ac- 
cretion disk in RUKL-type models at the canonical values given 
in RUKL05 and UKRL06. The RUKL propeller luminosity L is 
then obtained in a straightforward manner as 



L 3 5*53*f» 4 (P v6l /0.01sr 3 /4o 1 . 



(8) 



In Eq. ©, the units of all variables are as before. 

Comparison of Eqs. (01 and © immediately leads to the fol- 
lowing conclusions about IS -type and RUKL-type propeller lu- 
minosities. First, the scalings with p and M are almost identi- 
cal for the two types. Secondly, the scaling with the neutron-star 
spin period P S pm is stronger (-3 instead of -2) for the RUKL-type 
than for the IS-type. Finally, for identical values of the variables 
p, M, and P sp m, the RUKL-type propeller luminosity is about 
three orders of magnitude lower than the IS-type propeller lu- 
minosity. Conversely, at fixed values of fi and M, roughly equal 
luminosities are given by the two types if the spin-rate for the 
RUKL type is about an order of magnitude higher than that for 
the IS type. 



As indicated earlier, in this work we are exploring the prop- 
erties of such propellers as described above during the relatively 
early stages of post-supernova binaries containing pre-LMXBs, 
when the binary orbits are expected to be appreciably eccentric, 
as explained in Sec. 16.11 In such a system, the mass-transfer rate 
M, r through the inner Lagrangian point L\ is expected to vary pe- 
riodically with the orbital phase, as detailed below in Sec. [3] This 
flow of matter forms an accretion disk because of its large spe- 
cific angular momentum, as explained above, and slow viscous 
effects in the disk modify the profile of the above periodic mod- 
ulation (making it less sharp), and the resultant periodic profile 
is that which is shown by the mass-arrival rate M at the neutron 
star. The propeller luminosity then follows suite, showing a pe- 
riodic modulation, as described by Eq. ([3} for the IS-type torque 
or Eq. ([8]) for the RUKL-type torque. In this scenario, therefore, 
we identify the 6.67 hour period of IE with the binary period of 
a young, eccentric pre-LMXB, which is expected to turn much 
later into a standard LMXB after passing through further inter- 
mediate phases (see Sec. |6.41 i. In the next section, we give details 
of the expected nature of the mass-transfer modulation M tr {8) at 
the orbital period. 



3. Orbital modulation of mass transfer 

The problem of orbital modulation of mass transfer in eccen- 
tric orbits has been studied by a number of authors over almost 
three decades now, adopting various approaches appropriate for 
various aspects of the problem they have studied. These aspects 
have covered a considerable range, from a scrutiny of the con- 
cept of the Roche lobe in an eccentric orbit flAvni 1 976), to a 
study of test-particle motion through numerical integration of 
the restricted three-body problem at or near periastron passage 
dLubow & Shu 1975l l, to explicit calculations of orbital phase- 
dependent flow through L\ from a suitably-modeled stellar en- 
velope ( |Joss & Rappapo rt 1984| and references therein). For our 
purposes here, we have adopted the results of the calculations 
described by Brown and Boyle ( |Brown & Boyle 1984| hereafter 
BB): these authors described the flow through L\ from the at- 
mosphere of the lobe-filling companion with a scale height H as 
a sort of nozzle flow through the inner Lagrangian point, inte- 
grating over a Maxwellian distribution of velocities (character- 
ized by thermal velocity scale vj) for the stellar matter. Their 
final result for the rate of mass transfer as a function of the true 
anomaly 8 is given by: 



M tr {0) = M 



y 



1 +, 



y p 1 + e cos 8 



exp 



-yfje 



1 - cos 8 
1 + e cos 8 



(9) 



In Eq. (0, e is the orbital eccentricity, and the dimensionless 
function y{8) is the ratio of the phase-dependent equivalent 
Roche-lobe radius R(8) of the companion to the phase-dependent 
orbital distance d(8) in the eccentric orbit, y p being the value of 
y at periastron (8 = 0). From standard geometry of ellipses, d(8) 
is given by: 



d(8) 



1 +, 



1 + e cos ( 



(10) 



where p = a(l - e) is periastron distance. Finally, fi = p/H is 
the the dimensionless scale-height parameter introduced by BB. 

It is convenient to work in terms of the ratio y as it varies rel- 
atively slowly with orbital phase (and is, in fact, independent of 
this phase for a non-rotating companion: see below). The other 
properties of the binary system that y depends on are (a) the mass 
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ratio Q = M c /M x , M c being the mass of the low-mass compan- 
ion, and (b) the rate of rotation Q. c of the companion, usually ex- 
pressed in units of the orbital angular velocity Q p at periastron as 
A = Q. c IQ. p . The scale Mq = ^TjtjppHvjpQ of the mass-transfer 
rate in Eq. (0 is set by the above velocity scale v T , the scale-size 
pH for the effective cross-section of the above "nozzle", and the 
basic density scale po in the stellar atmosphere. 

Prescriptions for y have been given in the 1970s 
and '80s; we use here the generalized Joss-Rappaport 
( |Joss & Rap paport 1984) expressions adopted by BB, namely, 



y = A-Blogg + COog Qf, 

where the coefficients in y are given by: 

A = 0.398 - 0.026/: + 0.004ZT 3/2 ] 
B = -0.264 + 0.052/T - 0.015tf 3/2 \ 
C = -0.023 - 0.005/T J 



(11) 



(12) 



and the variable K depends on the above rotation parameter A 
and the orbital phase as: 



K = A l 



d+ef 



(1 + ecosdy 



From Eqs. ([TTll-dTSl. it is clear that, for a non-rotating com- 
panion with K = 0, y is independent of the orbital phase, and 
depends only on the mass ratio Q. Thus, for a given Q, R(ff) sim- 
ply scales with d(8) as the eccentric orbit is traversed. It is stellar 
rotation which modifies the Roche potential in such a way that 
this simple scaling is broken, and y depends on orbital phase. 
The phase-dependent factor in K goes back to the original work 
of Avni (1976). In our present work, we study the limits of (a) 
no stellar rotation, A — 0, and (b) synchronous stellar rotation, 
A — 1, to cover a range of possibilities (see below). The esti- 
mated accuracy in the above prescription for determining equiv- 
alent Roche-lobe radii is ~ 2%. 

Detailed models with the mass transfer profile M, r {6) given 
by Eq. (0 are described below. From general considerations, 
it is clear that this profile peaks at the periastron and that the 
sharpness of the peak depends on the quantity yBe. Since y ~ 1, 
and typical values of B for the current problem are in the range 
10 2 - 10 3 (BB), we see that the profile is expected to be sharply 
peaked at the periastron even for realtively low values of eccen- 
tricity, such as e ~ 0.2. 

4. Viscous flow in accretion disks 

Matter transferred through L\ into the Roche lobe of the neutron 
star first forms a ring around the neutron star, the radius r r i ng of 
this ring being related to the specific angular momentum ( tr of 
the transferred matter as (Prin gle 1981[ l: 



' ring 



(14) 



Through effective viscous stresses, this ring spreads into an 
accretion disk, wherein matter slowly spirals inward towards the 
neutron star as the viscous stresses remove angular momentum 
from it. The accretion disk extends from its outermost radius r out 
inward upto the magnetospheric boundary r m , where the pro- 
peller torques expel the matter by depositing energy and angular 
momentum in it, as explained above. 

The rate M at which the matter drifting radially inward 
through the accretion disk arrives at r m depends, therefore, both 
on the profile of mass supply M, r at L\ , as described above, and 



on the rate of viscous radial drift through the accretion disk, 
which occurs on a timescale f njc . 

In a quasi-steady state, the relation between the two profiles 
M(t) and M rr (t) is of the form 



M(t) 



f 

Jt-N*P„rb 



M tr (t )f(T)dt 



where 



to 



tvi si 



(15) 



The convolution integral in Eq. ( fT3l ) describes the viscous drift 
with the timescale t v i SC of the mass supplied to the disk at earlier 
times fo at the rate M, r (to), as indicated above. In principle, the 
integral extends over all previous times, but in practice it is suffi- 
cient to keep track of only about N orbital periods in the past (as 
the lower limit of integration indicates). This is so because of the 
rapid fall of of the viscous-evolution profile f(f) of the accretion 
disk at large values of r (see below). 

Viscous-evolution profiles have been calculated analytically 
and numerically at various levels of approximation by several 
authors ( |Lynden-Bell & P ringle 1974 Ligh tman 1974[ ). For our 
purposes here, we have adopted an analytic approximation of the 
generic form 



(13) f(?) = t 



(16) 



introduced and utilized by Pravdo and Ghosh 
(IPravdo & Ghosh 20011 hereafter PG). This reference has 
discussions of earlier analytical and numerical investigations. 
The generic PG profile in Eq. (fT6b reaches its maximum at 
t = 1, and decays subsequently as r~". Clearly, therefore, most 
of the contribution to the above convolution integral comes 
from those orbital cycles which are closest to the earlier time 
t o = t — tyisc, and N is determined by the sharpness of the fall 
of the profile, i.e., n. In our computations, we estimated the 
optimal values of N by running test cases with increasing values 
of N until the desired accuracy was obtained. For example, in 
the best-fit case reported below, we found that N - 9 gave an 
accuracy of w 10%, while N = 15 gave an accuracy of m 1%. 
Given the error bars on the data points in the observed light 
curve, further accuracy was unnecessary. 

The following generic feature of viscous evolution of accre- 
tion disks is a key aspect of the phenomenon we are exploring 
here. Whereas the orbital modulation of the mass-supply rate 
M, r (t) to the disk at its outer radius r out is expected to sharply 
peaked at periastron for typical values of the scale height in the 
companion's atmosphere, as above, the viscous drift of matter 
through the accretion disk would decrease the sharpness of this 
modulation, since variations on timescales much shorter than 
t V i SC tend to be "washed out" by viscous diffusion. This is what 
makes the orbital modulation of the mass-arrival rate M(t) at the 
disk's inner radius r m gentler, and therefore also the modulation 
of the propeller luminosity L(t), leading naturally to light curves 
of the form observed in IE. Quantitative details follow. 

5. Model light curves 

We constructed model light curves for IE by combining the 
model of mass transfer described in Sec. [3] with that of viscous 
flow through the accretion disk described in Sec. [4] We then 
fitted these models to the observed light curve of IE in 2005 
(dL06). The fitting parameters were (J3, e), which come from the 
above BB mass-transfer model in elliptic orbit, and also (? vi - JC , 
n), which come from the above PG parametrized description of 
viscous evolution of accretion disks. In this introductory work, 
we kept B constant at a canonical value of B = 10 3 (BB), and 
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Fig. 1. X-ray light curve of IE. Shown is the observed light curve from dL06, superposed on the (common) best-fit model light 
curve for IS-type and RUKL-type propellers. Left panel: model curve for A — (nonrotating companion). Right panel: same for 
A — 1 (synchronously rotating companion). 



varied the parameters e, f l(JC , and n to obtain acceptable fits. For 
the viscous timescale, we found it more convenient to work in 
terms of the ratio k = t V j SC /P or t of this timescale to the known 
period P or i, m 6.67 hr of the system, which we of course identify 
with the orbital period in this model. The ratio k is of immedi- 
ate physical significance, since it measures the relative impor- 
tance of viscous diffusion to orbital modulation in the system. 
For k <K 1, viscous diffusion would be so rapid as to enable 
the disk flow to adjust to the orbital modulation of mass-supply 
rate, and flow-rate would essentially follow the supply rate. For 
k » 1, on the other hand, the viscous diffusion would be so slow 
as to wash out any rapid variations in the mass-supply rate, and 
the modulation would be essentially determined by the disk vis- 
cosity. As we see below, k values of a few seem to describe the 
IE system, indicating comparable importance of the two effects 
in this system. 



Table 1. Best Fit Model Parameters: IS-type torque 



Parameter 


Best fit value (A=0) 


Best fit value (A=l) 


K 


2.6 


2.6 


viscous-profile index n 


5.04 


5.04 


eccentricity 


0.405 


0.406 


x 2 


1.013 


1.012 



Table 2. Best Fit Model Parameters: RUKL-type torque 



Parameter 


Best fit value (A=0) 


Best fit value (A=l) 


K 


2.6 


2.6 


viscous-profile index n 


5.02 


5.02 


eccentricity 


0.400 


0.401 


x 1 


1.003 


1.006 



We fitted model light curves corresponding to both IS-type 
and RUKL-type torques to the data on IE, the best-fit values 
of the parameters being given in Table Q] and Table [2] In each 
case, we have considered both a non-rotating secondary (A = 0) 
and a synchronously-rotating secondary (A = 1), as indicated. 
Note that the best-fit values for the two types of torques are 
very close to each other, as may have been expected. This is 
so because the closeness of the scaling of L with M between the 
two types, as discussed in Sec. [2] since only this aspect of the 



torque is relevant for fitting the profile of the light curve. Other 
aspects, e.g., the fact that the RUKL-type propeller luminosity is 
about three orders magnitude below the IS-type propeller lum- 
nosity for identical vaues of fi, M and a>, have important conse- 
quences elsewhere, as detailed in Sec. 16.31 but not in this matter. 
Further, the absolute values of the observed luminosities in the 
light curves are easily accounted for, e.g., by having the stel- 
lar spin rate higher for RUKL-type torques by about a factor 
of 10 than that for IS-type torques, and p, M identical for the 
two types, as the scalings in Eqs. © and (0 show. This implies 



neutron-star spin periods in the range P s 



0.01 - 0.1 s, i.e., 



the canonical range for propellers, for both type of torques, as 
explained in Sec. Q] 

This closeness of best-fit parameters is reflected in the best- 
fit light curves, which are visually essentially identical for the 
two types of torques. In Fig.Q] we display this common best-fit 
light curve, superposed on the data on IE. 

Our inferred best-fit value of k in the above tables indicates 
that the dominant contribution to the convolution integral de- 
scribed in the last section comes from the second and third or- 
bits preceding the time of observation. The corresponding vis- 
cous timescale ?,,,■„. w 17.3 hr is consistent with a rather thick 
disk with h/r ~ 0.1 - 0.5 and a canonical value ~ 0.1 - 1 for the 



disk viscosity parameter (Shakura & Sunyaev 1973 1. This seems 
consistent with the results of the RUKL numerical simulations. 
Note also that the best-fit value of the viscous-profile index n 
is consistent with the range of values n ~ 4 - 5 generally ex- 
pected for neutron-star systems, as per the discussion given in 
PG. Indeed, we found that values of n in the above range gener- 
ally worked for the IE system. Regarding the orbital eccentricity 
e, the best-fit values are as given in the tables, and we found that 
values of the eccentricity e in the range ~ 0.35 - 0.45 genearlly 
worked for the IE system: we discuss this in the next section. It is 
clear, therefore, that the model explored in this paper can account 
quantitatively for the observed IE light curve in 2005, for both 
IS-type and RUKL-type torques. We discuss in Sec. |7.2| possible 
reasons for the apparently different, "jagged" light curve hinted 
at by the 2001 observations of this system (dL06). 

6. Formation & evolution of prototype systems 

As indicated in Sec. Q] we are exploring in this work a model 
for systems like IE wherein the binary system of a He-star and 
a low-mass star (left after completion of the CE evolution phase 
in which the extensive envelope of the evolved primary has been 
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expelled and its He-core left behind) produces the pre-LMXB 
when the He-star explodes in a supernova (SN), leading to a 
newborn neutron star with a low-mass companion. Essential fea- 
tures of the formation and subsequent evolution of such systems 
are, therefore, essential components of this model. We now dis- 
cuss these features in brief, considering in this section first the 
immediate post-SN status of the system, and then the evolution 
of this system with the low-mass companion in an eccentric orbit 
at or near the point of Roche-lobe contact at periastron, produc- 
ing a system like IE where orbitally-modulated mass transfer 
proceeds through the inner Lagrangian point, and the newborn, 
fast-spinning neutron star is operating in the propeller regime, 
expelling this matter instead of accreting. Subsequently, we sum- 
marize further evolution of such systems. 

6.1. Immediate post-SN systems 

A major question that concerns us here is the expected ec- 
centricity of systems formed by the SN in the above scenario, 
since this eccentricity is crucial for the proposed mechanism. 
Qualitatively, it is obvious that the immediate post-SN system is 
almost guaranteed to be highly eccentric, as the mass loss from a 
typical pre-SN system of, say, a 3M He-star and a M c = O.4M 
low-mass companion (see below) in forming the post-SN sys- 
tem of M x = l.4M e neutron star with its M c = O.4M low- 
mass companion is 1 .6M , which is close enough to maximum 
allowed value of mass loss ( = half of the initial total mass of 
3.4M for zero kick velocity) to ensure that the post-SN orbit 
would be very eccentric. We shall use these values for the stellar 
masses throughout the rest of this paper. 

To see this quantitatively, we can adapt the extensive cal- 
culations of Kalogera, who computed the probability of the 
formation of X-ray binaries as a funtion of orbital parameters 
(Kalogera 1996). In the following, we shall use the same masses 
for the pre- and post-SN system as given above for illustrative 
purposes. The probability density from Kalogera's work is: 



G(a,e) = 



Hi 



2ne 



[a(l-e 2 )] 



2^1/2 



(17) 



1 + e \l - 



exp 



2£- 



1 (2a -I , 
f + 1 



-1/2 



Here, 



z = 



1/2 



and I„ is the modified Bessel function of zeroth order. Further, a 
is the ratio of semimajor axes of the pre- and post-SN orbits, £ 
is the ratio of the total mass of the post-SN binary to that of the 
pre-SN one, and £ = cr/V r , cr being the velocity dispersion in 
the SN kick-velocity, and V r the orbital velocity of the explod- 
ing star relative to its low-mass companion just before the SN 
(Kal ogera 1996) . 

In the problem we are studying here, the semimajor axis of 
the post-SN binary is determined by Kepler's third law from our 
assumed stellar masses above, and the known orbital period of 
IE. However, when there is a kick associated with the SN, the 
inferred semimajor axis of the pre-SN binary is not determined 
uniquely by the semimajor axis and the eccentricity of the post- 
SN binary: rather, there is a range of values corresponding to the 




0.4 0.5 0.6 

Eccentricity 



Fig. 2. Formation probability-density Q(e) of immediate post- 
SN binaries as a function of eccentricity e for various values 
of the dispersion cr in the SN kick velocity (see text). Curves la- 
beled by the value of vks = cr in units of 100 km s . Each curve 
so normalized that J Q(e)de = 1 . 



range of the kick-velocity. Thus, there is a range in the values of 
a: it is well-known that the allowed range for a is limited from 
1/(1 +e) to 1 /( 1 -e), these limits being first identified by Flannery 
and van den Heuvel (1975). Thus, for our purposes, it is aprro- 
priate to integrate G(a, e) over the above allowed range of a, and 
display the resultant probability density Q(e) = J G(a, e)da as a 
function of the eccentricity e. We show this in Fig.|2]for various 
typical values of cr as indicated. In this figure, we have used the 
symbol vks there to denote cr in units of 10 5 m s _1 = 100 km s _1 , 
the typical scale for the SN kick dispersion, and we have normal- 
ized the probability density Q(e) so that J Q(e)de — 1 in each 
case. As explained above, the closeness of the value of ( « 0.53 
in this typical case to its lower limit for no binary destruction 
in the SN (this limit is 0.5 for zero kick velocity) ensures that 
the probability peaks at a high value of e, as Fig. [2] shows. It is 
clear, therefore, that such a pre-LMXB would generically have a 
considerable eccentricity at the time of its formation in the SN. 



6.2. Tidal-evolution phase of pre-LMXBs 

The above newly-formed pre-LMXB undergoes tidal evolution, 
wherein three simultaneous processes occur, namely, (1) tidal 
circularization, i.e., decrease in the orbital eccentricity <?, (2) tidal 
orbit-shrinkage or hardening, i.e., decrease in the orbital semi- 
major axis a, and (3) tidal synchronization, whereby the rota- 
tion frequency Q e of the low-mass companion approaches the 
orbital angular frequency Q = 27r/P or /,. These processes hap- 
pen through tidal torques, and their quantitative descriptions pi- 
oneered by Zahn (1977, 1978) are widely used for calculations: 
we use them here, as have P08. Complete equations are given in 
Zahn (1977), and an Erratum was published by Zahn (1978). We 
have found a further algebraic or transcription error in the orig- 
inal paper, which we describe below, and which seems to have 
gone unnoticed so far. 

Complete formulations for the rates of change of e, a, and Q c 
are given in Zahn (1977), but for our work here we shall utilize 
a widely-used simplification which comes naturally out of these 
formulations, namely, that the timescale for tidal synchroniza- 
tion comes out to be much shorter than that for tidal circular- 
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ization and tidal hardening (see, e.g. JMeibom & Mathieu 2 005, 
P08). This is appropriate, since we shall be interested in this 
work only in phenomena which occur on the timescales of tidal 
circularizatuion or longer. Under such circumstances, we can 
look upon the system as being roughly synchronous at all times, 
and describe tidal circularization and tidal hardening respec- 
tively by Zahn's (1977) equation (4.7) and the appropriately 
simplified (i.e., synchronized) version of Zahn's equation (4.3), 
thereby obtaining: 



e at 4 tf \a ) 



and 

1 da 



Ida k 2 IR\ 

-— = 114-^(1 +q) - 

a dt tF \al 



(18) 



(19) 



In equations ( fT~8b and ( fT9l ), q = 1 /Q in terms of the mass ratio 
Q = M c /M x defined above in Sec[3] k 2 is the apsidal motion con- 
stant for the low-mass companion, and tF is the "friction time" 
of Zahn (1977), which, for stars with convective envelopes (as in 
the present case) is given by Zahn's (1977) pioneering prescrip- 
tion of the turbulent eddy-viscosity timescale tEV- 



t F ~ t EV = (M C R 2 /L C ) 



1/3 



(20) 



Equations ( fT8l and ( fT~9b describe simultaneous tidal circu- 
larization and hardening of close binaries, but before presenting 
our results, we need to correct two errors related to them. First, 
if we define a circularization timescale t cm = -ej(dejdt) in the 
usual way, we get from eq. ( fTSl the result: 



1 



■ I 



M r R 



2\l/3 



63 k 2 q(l + q) \ L c 



(!)*■ 



(21) 



which would be identical to Zahn's (1977) equation (4.13), ex- 
cept that the factor of 4 on the right-hand side is missing in Zahn 
(1977). Unfortunately, this error has propagated over the years 
into numerous papers, e.g., in P08, in their equation (2X\ We 
have corrected this now. Secondly, in an erratum published in 
1978, Zahn corrected a few other (generally smaller) numerical 
errors, of which the one relevant to our work is that the numeri- 
cal coefficient on the right-hand side of our eq. ( fT8l should be 21 
instead of 63/4. In all calculations reported here, we have made 
these corrections. 

We have integrated eqs. ( fT8l and ( fT9l numerically for close 
binary systems like IE, with values of initial post-SN semimajor 
axes and eccentricities, a -, and e,-, chosen over a range of plausi- 
ble values for such systems. We find that, in all cases, the sys- 
tems circularize and harden in a way that, in the (e vs. a) plane, 
the circularization point is approached in a "cut off" like manner. 
This is shown in Fig.|3]for a possible prototype IE-like system, 
so chosen that the parameters of it evolve to those roughly corre- 
sponding to IE in ~ 2000 years. This cut-off approach is similar 
to what Meibom and Mathieu (2005) found. Of course, our de- 
tailed shape is slightly different from that of these authors, since 
they fitted their results to an assumed parameterized distribution 
shape applicable to observations on a collection of "normal" bi- 
naries. These details wil be given in a separate publication. For 



1 Because of this, the parameters adopted for IE by P08 and by our- 
selves in this work actually give /„■„. ~ 10 4 yr. In our work, we have 
used the values of the apsidal-motion constant ki given by Landin et 
al. (2009). 





1.5 

semimajor axis 



Fig. 3. Tidal evolution of a prototype IE-like system in the e vs. 
a plane. Semimajor axis a in units of solar radius. Note the "cut 
off" like approach to the circularization point (see text). 

our purposes here, we note that the total time t ci > c taken to reach 
this circularization point (Meibom & Mathieu 2005) can be ex- 
pressed roughly as: 



T() 



-2.55 



(22) 



where a, and <?,■ are the initial semimajor axis and eccentricity of 
the immediate post-SN orbit, and the scale parameter to is given 
by: 



T() 



1 (M C R 2 C 
2*2?(l+?)\ U: 



1/3 



(23) 



Equation ( 1221 is a rough analytic fit to the mumerical results, 
adequate for our purposes. Note that the scale parameter to de- 
pends on the companion mass M c , its value being to ~ 1 yr for 
the inferred companion mass of IE. 

It is clear from eq. (1221 ) that circularization is faster for orbits 
which are born more compact and more eccentric. The scaling 
with a is straightforward from the above equations of tidal evolu- 
tion; the scaling with e is more complicated (although inspection 
of the same equations gives some clue), involving details of the 
numerical solution. 

The lifetime r C! > c of the eccentric phase of the pre-LMXB is 
obviously also the lifetime of its orbital-modulation phase which 
we are investigating in this work. The sensitive dependence of 
this lifetime on the initial post-SN orbital parameters and the 
companion mass (through the scale parameter to and due to the 
mass-dependence of R c in eq. l22l ) makes for a wide range of 
possible values of this lifetime, ~ 10 3 - 10 8 years. 

A crucial point is, of course, that if the companion is at or 
close to filling its Roche lobe at periastron in the post-SN orbit, it 
must remain so throughout most of this eccentric phase in order 
for the scenario to be self-consistent. The size of the Roche lobe 
at periastron is simply p — a(l - e) multiplied by a well-known 
function of the mass ratio q. Since the latter does not change 
significantly during this phase, we need only study the evolution 
of the former. Our integration of the tidal-evolution equations 
show that, while a and e both decrease during this phase, p - 
a(l — e) decreases slowly through most of this phase, reaching a 
minimum and increasing thereafter at late stages. This is shown 
in Fig.|4]for the prototype IE-like system displayed in Fig.[3](see 
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above). Thus, if the companion is initially at or close to filling 
its Roche lobe at periastron, it will remain so over most of this 
phase, and if it is inside its Roche lobe initially, it is likely to fill 
its Roche lobe later during this phase. It is also seen that Roche- 
lobe contact ends at the last parts of this phase (when the orbit is 
nearly circular), since p increases and becomes roughly constant 
there. 





3000 
time 



Fig. 4. Evolution of periastron distance p = a{ \ — e) during tidal 
evolution of a prototype IE-like system (see text). Shown is p in 
units of the solar radius vs. time in years. 



Thus, this tidal-evolution phase is a rough measure of the 
lifetime of Roche-lobe contact and orbital modulation of the pro- 
peller output. After this, the pre-LMXB becomes detached, and 
remains so until angular-momentum loss through gravitational 
radiation and/or magnetic braking brings it back to Roche-lobe 
contact on a long timescale of 10 s - 10 9 yrs. We discuss this 
phase below in Sec. 16.41 



6.3. Duration of propeller phase 

When the above tidal-evolution phase ends, is the neutron star 
still operating in the propeller phase? To answer this question, 
we consider the spindown of the neutron star from an initial spin 
period P\ pin to a final, longer spin period P{ pin under the action of 
the propeller torque given by either the IS-type torque (Eq. JT]) 
or the RUKL-type torque (Eq. Q). In each case, this spindown 
is decsribed by 



to N 1 



tL) ItD t, 



(24) 



prop 



where / is the moment of inertia of the neutron star and t prop is 
the propeller spindown timescale. 

First consider IS-type torques, for which t prop is given by 

r-(2GM ) 2/7 I c -3 _8 2 

t prop = 3 V2 ^ 8/7a T 3/7 * 2.5 x \Q s Mjp^mlh 5 yr, (25) 

where I45 is / in units of 10 45 gm cm 2 , and other units are as 
before. Equation ( f24b can be integrated readily in this case, the 



total spindown time T prop from P' spin to P\ pin being: 



r prop = 2.303 t P ro P lo E {P f SB JP\ pin ) 



As discussed earlier, the ratio P spi „/P' spin is believed to be in the 
range 10- 100 (Ghosh 1995 and references therein), and its exact 
value does not matter because of the logarithmic dependence. On 
taking mx = 1 .4 and the corresponding moment of inertia for a 
standard modern EOS, we arrive at 



2 x 10 6 yr 



(27) 



for canonical values of M, p. and /. 

Now consider RUKL-type torques, for which t pr0 p is given 

by 



tpnp = t (P S pi n /0.01s), where 
t * 2.3 x lO^'V^ 11 ^ yr. 



(28) 



Equation d24l i can be integrated readily in this case also, the total 
spindown time T prop from P\ pin to P S spin being: 

Tprop = to(P f spin ~ P spin )/(0.0U) - fo^LjO.Ols), (29) 



spin ' 



the second equality in the above equation coming from the fact 
that the ratio P[ pin IP l sp i n is believed to be in the range 10 - 100, 
as indicated above. The numerical value of T prop in this case is 
thus 



spin' spin' 



(26) 



Tp r op * 2.3 x lOXMo 1 W^ ; „/0.01s) yr, (30) 

which implies that, for canonical range P[ pin ~ 0.1 - 1 s, as 
indicated in Sec. [2] we arrive at 

T prop ~ 2 x 10 8 - 2 x 10 9 yr (31) 

for canonical values of the variables M, p and /. 

In comparing the total spindown times T pr0 p given by the two 
types of torques, we notice that the time taken by the RUKL- 
type torque is 2-3 orders of magnitude longer than that taken by 
the IS-type torque for identical values of M, p and /. This re- 
flects the relative weakness of the RUKL-type torque discussed 
in Sec. [2] Next, comparing the values of T prop given by the above 
two types of propeller torques with the lifetime r arc of the eccen- 
tric phase given in the previous section, we reach the following 
conclusions. For the IS-type torque, we find that, over most of 
the parameter space, the neutron star would still be in the pro- 
peller phase at the end of the above tidal-evolution phase of the 
binary. For the RUKL-type torque, we find that this conclusion 
is valid over the entire parameter space. Thus, the RUKL-type 
torque makes the conclusion stronger. 

As shown above, the companion has moved out of Roche- 
lobe contact by the time that the tidal-evolution phase of the bi- 
nary reaches conclusion, so that mass transfer stops, and so does 
the propeller action and its consequent soft X-ray production. 
Accordingly, throughout this first Roche-lobe contact phase, we 
expect the system to be in the propeller phase. 

6.4. Re-contact with Roche lobe & LMXB phase 

After orbit circularization and the loss of its first Roche-lobe 
contact, as described above, the pre-LMXB thus ceases to be 
an X-ray source. But its orbit shrinks (i.e., the binary hardens) 
on a long timescale (~ 10 s - 10 9 yr) due to two mechanisms of 
angular momentum loss from the system, viz. , graviational radia- 
tion and magnetic braking (Ghosh 2007 and references therein). 
These are the standard mechanisms through which short-period 
pre-LMXBs are believed to harden, until Roche-lobe contact is 
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regained and mass transfer restarts. But the transferred mass 
is now accreted by the neutron star, because its spin has been 
slowed down sufficiently over this long time that it acts as an ac- 
cretor and not a propeller at the (large) mass-transfer rates that 
occur at this second Roche-lobe contact in the circularized bi- 
nary. The system thus turns on as a canonical LMXB now, emit- 
ting strongly (L ~ 10 37 - 10 38 erg s _1 ) in the canonical X-ray 
band characteristic of emission from the neutron-star surface, 
rather than the soft X-ray band characteristic of propeller emis- 
sion from the vicinity of the magnetospheric boundary. 

The timescale ten of orbit shrinkage due to gravitational ra- 
diation is given by (see, e.g., Faulkner 1971, Banerjee & Ghosh 
2006): 



tcR ~ 2 x 10 9 



x m c \6 h j) yr 



(32) 



where mj = m x + m c , and all masses are in solar units. In this 
equation, we have scaled P or h to the value for IE, and substi- 
tution of the masses we have used above for this system gives 
tat ~ 4 x 10 9 yr. Generally, IE-like systems with shorter pe- 
riods and/or somewhat different companion masses will have 
tcR ~ 10 8 — 10 9 yr. Magnetic braking is believed to be compara- 
ble or weaker in strength to shrinkage by gravitational radiation 
at these orbital periods, so that the above estimate is a reasonable 
one for the IE-type systems we have in mind here. 

Thus, the system become a canonical, bright LMXB with a 
circular orbit and P or b in the range of, say, 2 - 10 hours. It is well- 
known that systems with P or b exceeding about 12 hours can- 
not come into Roche lobe contact by the above orbit-shrinkage 
mechanisms, since the time requied would exceed the Hubble 
time, as eq. ( 1321 readily shows. However, these long-period sys- 
tems also come into Roche-lobe contact eventually, as the low- 
mass companion completes its main-sequence evolution and ex- 
pands. These systems thus also become canonical long-period 
LMXBs with circular orbits. The lifetime of this standard, bright 
LMXB phase is t LMXB ~ 10 s - 10 9 yr. 



7. Discussion 

In this work, we have explored a pre-LMXB model of IE, 
wherein the eccentric orbit of the very young pre-LMXB causes 
an orbital modulation in the mass-transfer rate, and the new- 
born, fast-rotating neutron star operates in the propeller regime, 
the propeller emission in soft X-rays following the above 
modulation after viscous smoothening in the accretion disk. 
In this section, we first discuss first some essential spectral 
and luminosity-dependent features of IE, and their connec- 
tions with corresponding features in old, low-mass, soft X-ray 
transients (SXRTs) in their low/quiescent states, the prime ex- 



ample of this class being Aquila X-l ( |Campana et al. 1 998). 
Note that the well-known transient accretion-powered millisec- 
ond pulsar SAX J1808. 4-3658 also shows a similar behavior 
(Stell a et al. 20001 1. In these classes of low-mass X-ray binaries 
with old neutron stars, the neutron star is thought to operate in 
the propeller regime when the sources are in their low/quiescent 
states during decays of their outbursts. We then compare our 
model with the magnetar model which has been proposed re- 
cently for IE (P08), and discuss how distinction between the 
two kinds of models might be attempted in future. Finally, we 
summarize our conclusions. 



7.1. X-ray spectra 

The XMM-Newton/EPIC (0.5-8 keV) X-ray spectra of IE have 
been described by dL06. The time-averaged spectra from the 
2005 low-state observations, when the source luminosity was 
L ~ 10 33 erg s _1 , can be fitted by a two-component model con- 
sisting of a blackbody (BB) of temperature kTbb ~ 0.5 keV and 
an equivalent blackbody radius R/,b ~ 0.6 km, plus a power-law 
(PL) of index F ~ 3, with ~ 70% of the total flux coming from 
the blackbody component. Alternatively, the second component 
can also be a blackbody with a higher temperature. A re-analysis 
of the earlier 2001 XMM-Newton data, when IE had a higher 
luminosity (by a factor ~ 6), yielded a similar two-component 
(BB+PL) model with essentially the same blackbody tempera- 
ture kTbb and power-law index T, but a larger equivalent black- 
body radius Rbb ~ 1.3 km, and a higher contribution from the PL 
component (the blackbody contribution was ~ 50% of the total 
flux as opposed to the above ~ 70%), which made the overall 
spectrum harder (dL06). 

We stress the remarkable similarity of the above observations 
with those of the spectra of SXRTs in their low/quiescent states 
(when the neutron stars in them are believed to be functioning in 
the propeller regime), taking the well-known source Aquila X-l 
as the example. A detailed analysis of the BeppoSAX observa- 
tions of Aquila X-l in 1997 (Campan a et al. 1998] > has yielded 
the following results. At the lowest state, with source luminos- 
ity L ~ 0.6 x 10 33 erg s"\ the (BB+PL) fit had a BB of tem- 
perature kTbb ~ 0.3 keV and an equivalent blackbody radius 
Rbb ~ 0.8 km, plus a power-law (PL) of index r ~ 1, with ~ 60% 
of the total flux coming from the blackbody component. As the 
luminosity increased by a factor ~ 150 to L ~ 9 x 10 34 erg 
s _1 , the (BB+PL) fit yielded a BB of temperature kT h b ~ 0.4 
keV and an equivalent blackbody radius Rbb ~ 2.6 km, plus 
a power-law (PL) of index F ~ 1.9, with ~ 20% of the total 
flux coming from the blackbody component. Remembering that 
the total range of luminosities in these Aquila X-l low-state ob- 
servations during outburst decay is roughly 10 33 - 10 35 erg s _1 



(Campana et al. 1998| l, essentially identical to that of the IE ob- 
servations reported by dL06, the correspondence is very sugges- 
tive. 

SXRTs are believed to be old systems with a neutron star 
and a low-mass companion in a close circular orbit, under- 
going outbursts due to instabilties either in the accretion disk 
or in the mass supply from the low-mass companion. In their 
low/quiescent states during decays of outbursts, the fast-spining 
neutron star (spun up by accretion as per standard LMXB sce- 
nario) is believed to operate in the propeller regime. What we 
suggest in this work is that IE-like systems are very young sys- 
tems in the same regime: the young systems can show orbital 
modulation because of the orbital eccentricity, while the old sys- 
tems are in circular orbit and cannot show such orbital modula- 
tion. However, the spectral signatures are very similar at similar 
luminosities, which supports our basic suggestion. We note that 
the timescales associated with IE outburst appear to be ~ 2 - 3 
years while those associated with Aquila X-l outbursts appear 
to be ~ 30 - 70 days. It is possible that the basic phenomenon 
is rather similar in the two cases, and that the difference in detail 
is caused by the fact that accretion onto the neutron-star surface 
(with attendant high luminosities and hard X-ray spectra) does 
occur at the high states during the outbursts for old systems like 
Aquila X-l, but not for young systems like IE. 

A comprehensive theory of the emission spectra of pro- 
peller sources appears to be lacking, though Illarionov and 
co-authors have studied some effects of Comptonization 
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in propellers in wind-accreting 
pilarionov & Kompaneets 1990 



Attempts at constructing sue 



massive X-ray binaries 
llllarionov et al. 1993l >. 
theory for propellers in 



pre-LMXBs and in old LMXBs in low/quiescent states is clearly 
beyond the scope of this paper, and we shall confine ourselves 
here to the comment that the importance of Compton heating, 
considered in the above works on propellers in massive binaries, 
is also likely to be crucial for the systems we are focusing on in 
this work, as the observed power-law tails in the spectra at low 
luminosities suggest. These tails are particularly prominent in 
the low-state spectra of Aquila X-l (Campana et al. 1998). 



7.2. Luminosity dependence of light curve 

dL06 have compared the IE light curve in the 2005 low-state 
observations with that during the 2001 observations when the 
source luminosity was a factor ~ 6 higher. While the former 
light curve is relatively smooth with some cycle-to-cycle vari- 
ations, the latter one shows more complex, somewhat "jagged" 
structure, with an occasional dip. Further, the pulsed fraction de- 
creses from ~ 43% to ~ 12% as the luminosity increases. We 
discuss qualitatively how such features may arise. First, a pro- 
peller system is inherently more fluctuating than an accreting 
system, because of a variety of fluctuations possible at the site of 
shock-heating and outflow. As mass-supply rate through the ac- 
cretion disk increases, these fluctuations may increase, causing 
more complex profiles. Secondly, accretion disks in low-mass 
systems like LMXBs and CVs are thought to develop structures 
at their outer edges, which obscure emission from the compact 
object, and lead to dips. If these obscuring structures increase 
in size as mass-supply rate through the accretion disk increases, 
this would provide a natural explanation for the above appear- 
ance of the dips. Thirdly, as the mass-arrival rate M at r m in- 
creases, r m decreases (see Sec. 13, matter at the magnetospheric 
boundary becomes hotter, and the propeller becomes less super- 
sonic, ultimately becoming subsonic. Now, it is well-known that 
the subsonic propeller torque N su b ~ p?^l 2 JGM x is independent 
of M (see Ghosh 1995 and references therein), and so will not 
follow the modulations of M. Hence, as M and L increase, the 
following phenomenon is likely to happen. As the upper limit 
of the excursions in M goes beyond the critical cross-over point 
from supersonic to subsonic propeller regime, the pulsed frac- 
tion will decrease because that part of M which is above this 
critical point will not contribute to the pulsed flux, and this de- 
crease will increase with increasing M. This may be a natural 
explanation for the above observation of reduced pulsed fraction 
at higher luminosity. More quantitative considerations will be 
given elsewhere. 



7.3. Comparison with magnetar model 

In a recent paper, P08 have described a model in which IE is a 
magnetar, i.e., a neutron star with a superstrong magnetic field 
~ 10 15 G with a low-mass companion. The 6.7 hr period is in- 
terpreted in this model as the spin period of the neutron star, the 
idea being that a neutron star with such strong magnetic field 
as above can be spun down to such long spin period, or such 
low spin frequency, in ~ 2000 yrs. Magnetars are a fascinating 
possibility, and their relevance to soft gamma repeaters (SGRs) 
and possibly to anomalous X-ray pulsars (AXPs) has been the 
subject of much recent study. P08 have invoked an analogy with 
polars or AM Her-type cataclysmic variables containing white 
dwarfs with unusually strong magnetic fields, wherein torques 



acting on the magnetar spin it down in a short time to spin peri- 
ods in close synchronism with the binary orbital period. In this 
analogy, they have been inspired by the similarity of the shape 
the IE light curve to those of AM Her systems. 

We have desribed in this work a model which does not re- 
quire a neutron star with a superstrong magnetic field, but rather 
interprets the 6.7 hr period as the orbital period of the binary 
system consisting of a neutron star with a canonical magnetic 
field of ~ 10 12 G with a low-mass companion, the newborn, 
fast-rotating neutron star being in the propeller phase, and the 
propeller emission being modulated in the eccentric orbit of a 
young post-SN binary. We find that the observed IE light curve 
can be quantitatively accounted for by our model. Our analogy 
is with propeller regimes of SXRTs like Aquila X-l in their 
low/quiescent states, which we consider to be old, circularized 
analogues of IE which are no longer orbitally modulated, but 
which have remarkably similar spectral properties. In this anal- 
ogy, we have been inspired by the similarity between IE and the 
SXRTs in both the spectral characteristics and their changes with 
source luminosity, as well as the shapes of the outbursts and the 
way in which propeller-like properties emerge at low luminosi- 
ties during outburst decays. 

An interesting question is that of possible discriminators be- 
tween the above two models. It appears to us that if all observed 
properties of IE and similar systems can be accounted for by 
known characteristics of early stages of pre-LMXBs born ac- 
cording to the standard CE evolution and He-star supernova sce- 
nario, such as we have described in this paper (or by other pos- 
sible models involving standard evolutionary scenarios), there 
would not be any compelling need for invoking exotic objects 
like magnetars for this class of objects. On the other hand, if one 
finds unique observed features in this class of objects that cannot 
be explained at all within the framework of standard evolution- 
ary scenarios, presence of magnetars in such objects may well be 
hinted at. However, answering this question is beyond the scope 
of this paper: we are pursuing the matter, and the results will be 
reported elsewhere. 

7.4. Conclusions 

The work reported here suggests that IE-type systems are early 
stages of pre-LMXBs born in the SN of He-stars in binaries of 
(He-star + low-mass star) produced by common-envelope (CE) 
evolution. As long as the post-SN binary is eccentric, and the 
neutron star is in the propeller regime, soft X-ray emission mod- 
ulated at the orbital period may be expected to occur. As the orbit 
circularizes, modulation would stop, and as the low-mass com- 
panion moves out of Roche-lobe contact, the source would not 
be observed in X-rays. The companion would come into Roche- 
lobe contact again on a long timescale due to orbit shrinkage by 
emission of gravitational radiation and magnetic braking, and/or 
by the evolutionary expansion of the companion. This would 
lead to a standard LMXB: an old neutron star in circular orbit 
with a low-mass companion. Thus, steady-state arguments, with 
lifetimes of IE-type systems estimated at ~ 10 6 — 10 7 yrs and 
those of LMXBs estimated at ~ 10 s — 10 9 yrs, would lead us 
to expect ~ 1 IE-type systems per ~ 100 LMXBs, which is 
consistent with current observations. However, we must be care- 
ful here, as these are overall arguments for the whole popula- 
tion. If one specifically investigates young supernova remnants 
(SNRs), the chances of finding such systems may be consider- 
ably higher, since eccentric binary systems are to be found pref- 
erentially in such SNRs. More detailed considerations will be 
given elsewhere. 
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The lifetime of the eccentric -binary phase may be increased 
by an effect we have not included in this introductory work. The 
effect is that of an enhancement of eccentricity when mass and 
angular momentum are lost from a binary system which is al- 
ready eccentric. This dynamical effect is well-known in the lit- 
erature (see, e.g., Huang 1963 1 and its applications to compact 
X-ray binaries have been made earlier (Gh osh et al. 198 lT >. For 
an eccentric compact binary with the neutron star in the propeller 
regime leading to the loss of both mass and angular momentum 
from the system, such considerations are applicable. However, it 
is possible that, at the rates of mass transfer and loss inferred for 
IE-type systems, this effect is a minor one. 

Several lines of further investigation are naturally suggested 
by the considerations we have given in this paper. Foremost 
among them is a theory of the spectral characteristics of pro- 
peller emission in disk-fed propeller systems. This would help 
clarify the remarkable spectral similarity (including changes in 
spectral parameters with luminosity) between IE and SXRTs 
like Aquila X-l in their low/quiescent state, as described in 
Sec. 17.11 A search for point soft X-ray sources in other young 
SNRs would clarify the observational situation greatly. We note 
that these sources may or may not be periodically modulated, 
as we have argued in Sec. 17.21 that such modulations may de- 
crease and disappear in certain luminosity states. However, the 
spectral characteristics would still be a most valuable diagnostic. 
These and other investigations are under way, and results will be 
reported elsewhere. 
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